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FOREWORD 

The  research  reported  herein  was  conducted  in  the  Structural  Integrity 
Branch  of  the  Structural  Mechanics  Division,  for  the  Air  Force  Flight  Dyna- 
mics Laboratory,  Air  Force  Systems  Command,  Wright-Patterson  Air  Force  Base, 
Ohio,  under  project  1367,  "Structural  Integrity  for  Military  Aerospace  Vehi- 
cles," and  work  unit  13670336,  "Life  Analysis  and  Design  Methods  for  Aero- 
space Structures,"  by  Robert  L.  Neulieb.  The  research  was  conducted  from 
- September  1976  to  March  1977. 
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SECTION  I 
INTRODUCTION 

The  two-  or  three-parameter-Weibull  distribution  is  commonly  used  to 
describe  times  to  fatigue  failures.  Extensive  work  has  been  undertaken 
to  identify  the  shape  parameter  for  several  aircraft  materials  (References 
1 and  2).  Using  these  estimates  for  the  shape  parameter,  it  is  necessary 
to  obtain  the  scale  parameter.  Impellizzeri  et  al  (Reference  3)  and 
Freudenthal  (Reference  4)  have  shown  that  it  is  possible  to  obtain  in 
closed  form  a distribution  for  the  maximum  likelihood  estimate  of  the 
scale  parameter  divided  by  the  time  to  failure.  Freudenthal  (Reference 
4)  observed  that  in  terms  of  estimating  the  time  to  first  failure  in  a 
given  size  fleet,  little  change  in  reliability  is  obtained  for  the  same 
scatter  factor  by  conducting  more  than  one  test  to  estimate  the  scale 
parameter.  This  result  is  extremely  important  and  can  be  extended  to  the 
task  of  estimating  any  number  of  failures  as  long  as  the  quantity  time  di- 
vided by  the  scale  parameter  estimate,  all  raised  to  a power  given  by  the 
shape  parameter,  is  small.  This  is  the  most  important  region  in  terms  of 
aircraft  structural  integrity,  since  only  a few  failures  can  be  tolerated. 

However,  care  must  be  taken  in  interpreting  the  above  result.  Although 
the  probability  of  failure  (down  to  a given  ratio  of  the  scale  parameter  es- 
timate to  time)  is  little  changed  by  conducting  additional  tests,  the  value 
of  the  scale  parameter  estimate  may  vary  as  the  precision  of  the  estimate 
increases  with  the  number  of  tests.  Hence,  the  precision  and  value  of  the 
estimate  of  a time  to  first  or  subsequent  failure  (for  a constant  ratio) 
may  also  vary  with  the  data  base. 

In  this  report,  for  the  case  when  the  shape  parameter  is  known,  the 
probability  distribution  functions  and  asymptotic  functions  for  times  to 
failure  as  a ratio  of  the  maximum  likelihood  estimate  of  the  scale  para- 
meter are  derived.  The  influence  on  the  probability  of  failing  due  to  the 
size  of  the  data  base  is  investigated. 
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SECTION  II 

DISTRIBUTION  FUNCTION  FOR  TIME  TO  FAILURE 

In  both  Reference  3 and  4 derivations  were  obtained  of  distribution 
functions  for  scatter  factors  when  the  Weibull  shape  parameter  is  known. 

By  following  the  same  procedure,  it  is  possible  to  derive  the  distribution 
function  for  time  to  failure  of  an  item  when  given  the  maximum  likelihood 
estimate  of  the  scale  parameter.  The  times  to  failure,  given  the  exact  value 
of  the  scale  parameter,  will  be  assumed  to  be  two-parameter-Weibull  distri- 
buted. That  is,  it  is  assumed  that  times  to  failure  are  exactly  two-para- 
meter-Weibull  distributed;  however,  the  scale  parameter  is  not  known  because 
it  is  estimated  from  a finite  set  of  sample  data.  In  using  the  estimate  of 
the  scale  parameter,  the  distribution  of  times  to  failure  is  not  Weibull  and 
must  be  derived. 

The  cumulative  distribution  function  for  time  to  failure  will  be  de- 
rived by  multiplying  the  density  function  for  the  scale  parameter  and  the 
Weibull  density  function  for  time  to  failure.  This  product  is  the  joint  den- 
sity function  for  time  to  failure  and  the  scale  parameter  since  the  tests 
for  the  scale  parameter  are  independent  of  inservice  failures.  The  joint  den- 
sity function  will  be  integrated  over  the  entire  possible  range  of  the  scale 
parameter  to  obtain  the  density  function  for  time  to  failure  alone.  Then 
this  density  function  will  be  integrated  from  the  smallest  possible  time 
(zero)  to  the  time  of  interest  to  form  the  cumulative  distribution  function. 

Since  the  density  function  for  a function  of  the  scale  parameter  is 
known,  the  density  function  for  the  scale  parameter  will  be  obtained  by 
transforming  the  random  variable.  The  density  function  for  a transformed 
random  variable  is  related  to  the  original  density  function  by  the  expres- 
sion given  in  Reference  5.  This  expression  is 


fx(x)  = fw(w(x)) 


d w 
dx 


(D 
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The  probability  density  function  for  the  two-parameter-Weibull  distri- 
bution is 


f(y)  s 


(2) 


The  maximum  likelihood  point  estimator,  j3 , of  the  scale  parameter )S 
from  n tests  for  times,  y,  to  failure  is 


(3) 


where  a is  the  known  shape  parameter.  It  has  been  shown  (Reference  6)  that 
2n  [p/fi)  has  a chi-square  distribution  with  2n  degrees  of  freedom. 
Therefore  the  probability  density  function  is  given  by 


(4) 


where  T (n)  is  a complete  Gamma  Function. 


The  probability  density  function  for  the  random  variable,  /3 , can 
be  obtained  by  a transformation  using  Equation  1.  Hence, 


P 

we  can  obtain 


letting 


. (£f 


¥1 


and  observing  that 


dw 


-( 


f-  $ 

) j?' 


f(«- 


^r{n)  \$J 


(5) 
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This  is  the  probability  density  function  for  the  scale  parameter  given 

A A 

p and  n.  Hence,  if  p is  determined  from  n tests,  the  joint  probability 
density  function  for  time  to  failure  and  0 is 


0 t 


(6) 


A 

since  the  determination  of  (3  and  the  failure  of  a given  item  are  indepen- 
dent. The  density  function  for  times  to  failure  is  given  by 


f (y)  • J f (y, >3)  d/3 

O 


(7) 


or 


a2  ya-,nn  Aan  i \° 

,t,)  * -m- 13  m) 


an+a+l 


_ Aa 

(ya  + nfi  \ 

e'V  Ra  / 


£a  7 d/0  (8) 


-a  / \a+l 

Letting  u = 0~  and  hence  -ay-jz)  dj3sdu 


we  obtain 


i (n)  J 


(9) 


Letting  v = u (ya  + n f3a)  and  hence  dv  = (ya  + n)0a)  du  , we  find 


f (y)  « 


aya_l  nn  0 


Aan 


A 


r ( n ) (ya  + n/3  ) 


f vn 

a*n+i  J 


e_v  dv 


(10) 
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But 


j v"  e_v  dv  = n T( n ) 


so  that 


f (y)  * 


gnntl 

(y«+n  £T 


(11) 


(12) 


It  is  convenient  to  convert  this  density  function  to  a cumulative  distri- 
bution function. 


F 


(13) 


F (y ) = ann+l  pf 


Aan  -y  -a-' 


° (qa  + n£  ) 


dq 


(14) 


Letting  z = qa  + n/3~  and  hence  dz  = aqa_l  dq  , we  obtain 

y°+n)8a 

Artn  m ^ 

(15) 


ya+n)8a 

<»>  = "n*'£a7* 

nP 


( y/p ) +n 


(16) 
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This  function  is  the  probability  distribution  function  for  times  to 
failure  when  the  population  of  times  to  failure  is  Weibull,  the  shape  para- 
meter, a , is  known,  but  a finite  data  sample  of  size,  n,  is  used  to  obtain 

A 

the  scale  parameter  estimate,  p.  That  is,  this  function  is  the  probability 
that  an  item  will  fail  between  time  zero  and  time,  y,  under  the  conditions 
given  above.  This  distribution  function  can  be  used  to  determine  the 
probabilities  of  the  first  or  subsequent  failures  occurring  before  time  y. 
Methods  of  determining  the  probability  of  the  first  and  subsequent  failures 
given  the  probability  of  an  item  failing  in  some  time  increment  are  given  in 
Reference  7. 

An  expression  for  the  distribution  of  times  to  first  failure  in  "m" 
items  can  be  developed  directly.  The  probability  of  zero  failures  in  "m" 
independent  items  experiencing  identical  loadings  is  given  by  the  proba- 
bility that  one  item  does  not  fail  when  raised  to  the  power  given  by  the 
number  of  items,  m.  The  probability  of  an  item  not  failing  before  time,  y, 
is  one  minus  the  probability  that  the  item  fails.  Hence, 


p0<y> 


( q 

1 (y/£)a  + n 


(17) 


where  P (y)  is  the  probability  of  zero  failures  in  "m"  items  up  to  time, 
y,  or  equivalently,  the  probability  that  the  first  failure  will  occur  after 
time,  y.  Equation  17  is  comparable  to  but  not  identical  to  the  reliability 
function  developed  by  Freudenthal  in  Reference  4.  The  difference  arises  in 
Reference  4 from  the  use  of  a function  of  time,  the  scatter  factor,  instead 
of  time  itself  in  the  derivation  and  application.  Classically,  the  defini- 
tion of  reliability  involves  time  directly,  as  in  Equation  17.  See  Refer- 
ence 8 for  a more  complete  discussion.  Differences  between  the  two  formu- 
lations are  not  expected  to  be  significant,  at  least  for  aircraft  appli- 
cations. 
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SECTION  III 

DISCUSSION  OF  THE  DISTRIBUTION  FUNCTION 

The  cumulative  distribution  function  F(y)  in  Equation  16  gives  the 
probability  that  the  fatigue  life  of  a specific  item  is  less  than  or  equal 

A 

to  y when  p is  the  maximum  likelihood  estimate  of  the  scale  parameter  ob- 
tained from  n tests  and  a is  the  known  value  of  the  shape  parameter.  The 
failures  are  assumed  to  be  two-parameter-Weibull  distributed. 

The  distribution  differs  from  the  Weibull  distribution  because  the 
Wei  bull  scale  parameter  is  not  known  exactly.  As  the  number  of  tests  con- 
ducted to  determine  increases,  this  distribution  should  approach  the 
Weibull  distribution. 

By  rewriting  F(y)  = 1 - G(y),  the  limit  of  G(y)  as  n— •»  can  be 

obtained. 


Lim 
n— © 


G (y)  * Lim  ( 7^77 ) 

n—  V(y/£)%n  ' 

= Lim  ( 
n— ® \ 


i/n(y/)9  )°  + 1 


f 


Lim  e-n^n(,/n(y/^)a+l) 
n-»« 


Lim 

n-»® 


^n(i/n(y/)8)  +1 
e"  n"' 


(18) 
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The  argument  is  in  an  indeterminate  form-§-as  n— , hence  L'  Hospital's 
rule  may  be  applied.  Let  the  argument  equal  A.  Then 


Lim  A = -Lim 
n— n— « 


$-{£n  i/n(y/^)a+l  ) 


dn 


■aSr(--> 


n— >® 


/n(y/^)  + 


-n 


-2 


09) 


Hence  Lim  6 (y)  * e 

n— »® 


<y  t$f 


and 


F (y ) = I - e 


(y  'ha 


Thus,  as  more  information  about  the  scale  parameter  is  obtained,  this 

A o 

distribution  does  approach  the  Weibull  as  expected,  and  p goes  top  . 


The  region  where  (y/>S )a  «1  is  important  to  analyze,  because  this 
implies  that  there  is  a small  probability  of  any  given  item  failing  and 
hence,  that  only  a few  failures  occur  in  a fleet.  Since  tests  on  either 
full  scale  or  large  scale  items  are  very  expensive,  the  information  obtain- 
able from  one  test  should  be  analyzed.  Thus,  let  F(y)  for  n = 1 be  F^(y). 


F,  (y) 


t 

(y/£)  +i 


(20) 
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t 


* 


If 


/ A \a 

( y/£ ) « i 


F,  ( y ) « ( y/p  ) 


(21) 


Now  condsider  the  maximum  information  which  can  be  obtained  from 

A 

conducting  tests  to  determine  p . 


Lim 
n— ® 


F (y ) « I -e 


/ A\a 

-(y/0) 


(22) 


(W  V 

y/p  ) , we  find  that 

. a xa 

for  small  values  of  (y/p  ) » 


Lim 
n— *•« 


F(y)  w (y/P  ) 


(23) 


For  application  to  aircraft  failures, (y/p)  is  almost  always  small. 

The  asymptotic  distribution  seems  to  be  the  same,  independent  of  the 
number  of  tests  conducted  to  estimate  the  scale  parameter.  This  is  consis- 
tent with  the  observations  made  by  Freudenthal  (Reference  4).  However,  the 
quantity  need  not  be  identical  for  different  numbers  of  tests  used  to 
estimate  the  scale  parameter.  Hence,  for  a given  probability,  y need  not 
be  identical.  In  addition  to  Equation  5,  the  distribution  of  the  statistic 
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F (y) 
n'-7 ' 

~Hy) 


where  F (y)  is  given  by  Equation  16  and  F(y)  is  Weibull  can  indicate 


A 


the  significance  of  possible  variations  in  p . 

Fn(y) 

Let 


implying  that 


x = 


-(»//3) 


l-e 


P = y 


(24) 


(25) 


as  F (y)  is  independent  of  both  yS  and  x.  By  using  the  transformation  of  a 
random  variable.  Equation  1,  the  probability  density  function  for  the  ratio 
X can  be  obtained  from  Equation  5.  Observing  that 


t ■*(-*(■- W,“"(-1^)' 


Fn<y) 

v2 


(26) 


we  find  that 


r ( x ) = 


(If  V 


,an 


Fnjy) 
T(  n ) 


(27) 


This  is  the  probability  density  function  for  the  ratio  X.  The  values  of 
X may  range  from  F (y)  to  ® as  the  denominator  ranges  from  1 to  0.  This 
probability  density  function  can  be  integrated  to  a cumulative  distribution 
function 
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R (x)  = 


. £ .a 

/ Fn(y Kn\p^)  _l 

(4r  «) 


>n<y> 


Let 


0-*n- 


and  dg  8 dg 


-■(4)aof^r^(— r^r- 


Fn(y) 


6 »a 


(29) 


This  form  can  readily  be  integrated  for  n = 1 and  yields 


R,  (x)  = ( I — F,  (y)) 


/ ft  \a 

(p'  y) 


(30) 


A .a 


R , ( x ) « e“l/x  when  (y/£)  «l 


(31) 


Additional  R (x)'s  can  be  determined  by  integrating  by  parts,  or  by  writing 
n 

a general  expression  involving  the  ratio  of  incomplete  and  complete  Gamma 
Functions  (complement  of  the  Gamma  distribution). 

Let  g = e~p  and  dg  = - e-p  d p . Then 
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A .a 


R‘«'-(#r^4n(1-^)p""e'np(^rdp  (32) 


/ & »a 

Let  2 s pn  (p/y) 


and  dz 


= n(^/y) 


dp 


Then 


R (x ) = 


T'MUfa  )/„(.-^) 


zn~'  e_z  dz 


(33) 


R (x)  « 


r(n,-(/3/y)ai’n  (l-  ^T-)) 


T{  n ) 


(34) 


where  T (n,  x)  is  an  incomplete  Gamma  Function  as  given  in  Reference  9. 
For  n = 2 


R3(„  = (1-a(^„(1-^l))(l-^lj 


/ A \a 

(0/y) 


(35) 


See  Figure  1 for  plots  of  R-j(x),  Rg(x)  and  R®  (x).  Table  1 gives  the 
probability  that  X will  be  in  specified  ranges.  There  is  significant  varia- 
tion in  X about  the  value  1 for  n = 1.  However,  the  value  of  X is  more  fre- 
quently greater  than  1.  This  implies  that  in  these  cases  the  probability 
of  failure  is  overestimated.  As  n increases,  this  variability  in  X is  re- 
duced. Conducting  additional  tests  to  estimate  the  scale  parameter  does  pro- 
vide increased  information.  However,  the  value  of  this  information  must  be 
weighed  against  the  high  cost  of  obtaining  it. 
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(y/0)  = .01 
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TABLE  I 

PROBABILITIES  OF  Fp(y)/F(y)  BEING  IN  SPECIFIED  RANGES 


Range 

n = 1 

n = 2 

R(1 .25)  - R( ,8) 

.164 

.239 

R(2)  - R(.5) 

.473 

.647 

R(3)  - R(l/3) 

.669 

.840 

R(»)  - R(D 

.630 

.593 
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The  results  can  be  extended  to  the  case  of  a three-parameter-Weibull 
distribution  when  both  the  minimum  life  parameter  and  the  shape  parameter 
are  known.  When  the  minimum  life  parameter  is  known,  a parameter,  z,  the 
difference  between  times  to  failure  and  the  minimum  life  parameter  is  two- 
parameter-Weibull  distributed.  However,  in  terms  of  times  to  failure  this 
is  the  three-parameter-Weibull  distribution.  Thus  if  /jl  is  the  minimum  life 

z = y-/i  (36) 

F (z)  = l-e'(z/£)  (37) 

and 

-W) 

F (y ) = I - e (38) 


Hence,  when  /a  is  known,  a random  variable  which  is  three-parameter  distri- 
buted can  be  transformed  to  a variable  which  is  two-parameter-Weibull  distri- 
buted. All  properties  of  the  two-parameter-Weibull  distribution  thus  apply 
to  this  transformed  random  variable. 
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SECTION  IV 

CONCLUSIONS  AND  RECOMMENDATIONS 

The  probability  of  failure  in  the  region  of  small  probabilities  has 
the  same  form  whether  one  or  an  infinite  number  of  tests  is  used  to  esti- 
mate the  scale  parameter  when  the  shape  parameter  is  known.  This  con- 
firms Freudenthal ' s observation  (Reference  4)  that,  for  the  first  failure, 
scatter  factors  change  little  with  the  number  of  tests  conducted. 

More  information  can  be  obtained,  however,  by  conducting  additional 
tests  because  the  estimate  of  the  scale  parameter  can  vary  with  the  num- 
ber of  tests  conducted.  Equation  29  along  with  the  scatter  factor  or 
probability  of  failure  (Equation  16)  should  be  analyzed  to  determine  the 
number  of  tests  required. 

When  estimating  probabilities  of  first  or  early  failures  from  one 
or  more  tests  when  the  shape  parameter  is  believed  known.  Equation  16 
or  Freudenthal 1 s results  (Reference  4)  should  be  used. 
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